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O ABSTRACT 
Q 

^ Goodness-of-fit statistics are used to quantitatively establish the compatibility of CMB 

^ anisotropy predictions in a wide range of DMR-normalized, open and spatially-flat A, CDM 

QQ cosmogonies with the set of all presently available small-scale CMB anisotropy detection 

^ data. Conclusions regarding model viability depend sensitively on the prescription used 

^ to account for the la uncertainty in the assumed value of the DMR normalization, except 

for low-density, Oq ~ 0.3 - 0.4, open models which are compatible with the data for 
all prescriptions used. While large baryon-density (O^ ^ 0.0175/z~^), old (to ^15-16 

Q-i 

Q Gyr), low-density (Oq ~ 0.2 - 0.4), flat- A models might be incompatible, no model is 

incompatible with the data for all prescriptions. In fact, some open models seem to flt the 



data better than should be expected, and this might be an indication that some error bars 
are mildly overconservative. 

Subject headings: cosmic microwave background — cosmology: observations — large-scale 
structure of the universe 
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1. INTRODUCTION 

Recent and near-future measurements of CMB anisotropy, on a variety of angular 
scales, when used in conjunction with the predicted anisotropy in cosmogonical models, 
are in the process of transforming the CMB anisotropy field from that in which one tries 
to draw qualitative conclusions about the viability of broad-brush cosmological models, 
to that in which it will soon be possible to set quantitative constraints on parameters of 
some specific models and rule out other models. Until now all quantitative comparisons 
between model predictions and the data have made use of one of two simplifications: (1) 
data from one (or a few) experiments has been compared to predictions for one (or a 
class of) model(s); or, (2) data from a larger number of experiments has been compared 
to predictions for a single model. While clearly a necessary first step, this approach has 
led to a number of vague claims about the (in)compatibility of some model(s) with some 
subset of the data, which, while perhaps correct, need to be put on firmer ground. 

This work is a first attempt to compare all presently available CMB anisotropy de- 
tection data to predictions in a wide variety of observationally motivated cosmogonies, 
with the ultimate goal of deciding, in a quantitative manner, whether (or if) any of these 
models are compatible with the wealth of CMB anisotropy data. Such a quantitative ap- 
proach, using all available data, is essential if one wishes to draw robust conclusions about 
model viability. It will become more effective as the analyses are understood better, and 
as the data improves. The only other alternative is to wait a decade or so for a new CMB 
anisotropy satellite to address this issue. 

To qualitatively assess compatibility, Ratra et al. (1995, hereafter RBGS) and Ra- 
tra & Sugiyama (1995, hereafter RS) compared anisotropy predictions in 2 year DMR- 
normalized, gaussian, adiabatic, open and spatially-flat A, CDM models (with the values 
of Oo, and Qb chosen to satisfy non-CMB observational constraints, except in the 
flducial CDM case) to small-scale anisotropy data"^. Here we use these predictions, in 

It is important to bear in mind that a variety of statistical techniques and prescrip- 
tions (as well as different assumed CMB anisotropy spectra) have been used to determine 
the observational results (RBGS), that the usual prescription for accounting for calibration 
uncertainty, by adding it in quadrature, is not quite correct (RBGS), and that non-CMB 
contamination and/or subtraction might be an issue in some cases. 
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combination with a variety of goodness-of-fit statistics^, to quantitatively assess the com- 
patibihty of CMB anisotropy detections. In contrast to RBGS and RS, we explore more 
options for accounting for the Icr uncertainty in the DMR normalization, but in this pre- 
liminary analysis we ignore small-scale CMB anisotropy upper limits as well as the small 
correlations between data points from experiments with multiple windows (see §3). 

In a related analysis, Scott, Silk, & White (1995, hereafter SSW) used a Lorentzian 
approximation for the shape of the fiducial CDM model CMB anisotropy spectrum and 
concluded that it provided an adequate description of the anisotropy data (they took the 
data error bars to be symmetric). Here we use significantly more observational data, as well 
as revised estimates of some of the older data, and also use numerically computed CMB 
anisotropy spectra for a wider variety of models motivated by non-CMB observations. 
Consistent with the conclusion of SSW, for all prescriptions we have used to account for 
the allowed Icr range of the DMR normalization, low-density open models with Oq ~ 0.3 
- 0.4 are compatible with the data, so a CMB anisotropy spectrum that mildly rises to 
multipole moments / ~ 200 is compatible with the data^. 

2. SUMMARY OF COMPUTATION 

We consider 32 smaller-scale CMB anisotropy detections (almost entirely sensitive 
only to / < 200): FIRS (Ganga et al. 1994); Tenerife (Hancock et al. 1995); SK93 and 
individual-chop SK94 Ka and Q (Netterfield et al. 1996); SP94 Ka and Q (Gundersen et 
al. 1995); Python-G, -L, and -S (e.g., Piatt et al. 1995); ARGO (de Bernardis et al. 1994); 
MAX3, individual-channel MAX4, and MAX5 (e.g., Tanaka et al. 1995); MSAM92 and 
MSAM94 (Cheng et al. 1995); and WDHl (Griffin et al. 1995). For the i^^ detection, the 
observed bandtemperature STI^^~^ and the Icr upper and lower limits, <^2^(e)M *^"^(e)?' 
the DMR-normalized, 26 open and flat-A, model predictions for the bandtemperature Icr 
(gaussian) upper and lower limits, 5T^_^~^^ and 5T^_^~^j (range accounts for statistical and 

^ We use various such statistics since most observational error bars are asymmetric 
(i.e., nongaussian). 

^ As noted by SSW, such a rising CMB anisotropy spectrum is consistent with that ex- 
pected from radiation-pressure-induced oscillations at early times in the adiabatic structure 
formation picture. It is also consistent with that expected in versions of the isocurvature 
scenario. And there almost certainly are other models that are consistent with the data. 
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systematic uncertainty in the DMR normahzation, Stompor, Gorski, & Banday 1995), are 
given in RBGS (open models, DMR-galactic-frame normalization) and RS (flat-A models, 
DMR-ecliptic-frame normalization). Details may be found in these papers, and the model 
parameter values are listed in Table 1. 

To assess the effect of varying the DMR normalization, we consider 3 sets of model 
predictions 5T^_^~^^^ normalized to the lower Icr (= 5T^_^~^j)^ central (= [5T^_^~^^ + 5T^_^^^j]/2)^ 
and upper Icr (= 5TI^_^~^^) values of the DMR normalization. For each of these predictions 
we also either account for or ignore the model (DMR) "errors" cr^^_^~^ = [5T^_^~^^ — 5T^_^^^j]/2^ 
yielding 6 sets of predictions. Five prescriptions (/) are used to construct the corresponding 
observational numbers: (1) 5T,^ x = 5T,^ x and crj x = 5T,^ x — 5T,^ xi (2) 5T,^ x = 5T,^ x and 

V / (ejc (e) (e) (ejM (e) ' V / (ejc (e) 

^(e) = ['^^e). - <^^e)J/2; (3) ST^^^^ = [ST^^^^ + STl^^,]/2 and a^^ = [ST^^^^ - STl^^,]/2; (4) 
'^^e)c = <^^e) ^nd either = [ST^^^^ - ST^^^^] if ST^^^^ > ST^^^^^ or = [ST^^^^^ - ST^^^,] 
if ^Tl^)c < ^Tl)c^ and, (5) 5Tl^^^ = STl^^ and a^^ = STl^^ - STl^^, For each model, 
DMR-normalization value, model "error" prescription, detection and observational data 
prescription /, D'^j^ = [^Tl^e^c ~ <^^(m)c]/^* [where cr* is either = Icr^^^l or = ([('"(e)]^ + 
['^(111)]^)°^] i^ used as a measure of the deviation of the prediction from the observation. 
The -D^j-j are used to compute the reduced goodness-of-fit statistic xfj-) = S^=i [-^(j)]^/32 
(usual reduced for points drawn from a gaussian distribution). These are shown in Figs. 
1 - 6. 

3. DISCUSSION 

Some of the 32 detections we use here are not completely independent. As a result, 
there are slightly less than the 32 degrees of freedom we have assumed here, which, by 
itself, causes our reduced xfjj values to be slightly smaller than they should be. (A more 
accurate computation will require the appropriate correlation matrices.) One might hope 
to roughly compensate for this by focussing on the xfjj computed using -D^j-j with 
in the denominator [instead of ([c^*g-|]^ -|- [c^*jj^-)]^)°'^ which leads to an overestimate of the 
uncertainty due to cosmic variance (which has already been accounted for in the small-scale 
data error bars, and is an issue especially for data points from larger-angle experiments), 
and that due to systematic shifts in DMR normalization (which has already been accounted 
for by our use of 3 different DMR-normalization values)], but this would ignore the DMR 
noise uncertainty. 
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The skewness and kurtosis of the 780 -D^j-^ distributions, for each model, DMR- 
normalization value, observational data prescription, and model "error" prescription, is 
consistent with the range set by the variances of the skewness and kurtosis for 32 degrees 
of freedom drawn from a gaussian distribution. This means that less-compatible models 
(large xfjj in Figures) are less-compatible because of many somewhat deviant predictions, 
and not because of just a few extremely deviant predictions. 

Focussing on the nominal-DMR-normalized xfjj (Figs. 1 & 2), independent of the 
model "error" prescription, the only low-density (Oq ~ 0.2 - 0.4) flat-A models compatible 
with the data (i.e., with xfjj < 1.46, which is 2a [4.55% probability of xfjj being this large 
or larger] for 32 degrees of freedom drawn from a gaussian distribution) are the younger 
(to ~ 13 Gyr), lower baryon-density (O^ ^ 0.0075/z~^) ones^, while all open models are 
compatible with the data (in agreement with the qualitative conclusions of RBGS and 
RS). For the upper la value of the DMR normalization (Figs. 3 & 4), independent of 
the model "error" prescription, all flat-A models are incompatible, while low-density open 
models are compatible. At the lower la value of the DMR normalization (Figs. 5 & 6), 
the Oo = 0.1 open model is incompatible, most other models are compatible, and, in this 
case, conclusions regarding the viability of low-density (Oq ~ 0.2 - 0.4), old (to ~ 15 - 
16 Gyr), high baryon-density [Qb ~ 0.0175/z~^) flat-A models depend sensitively on the 
model "error" prescription^. 

Independent of the DMR-normalization value and model "error" prescription, low- 
density, Oq ~ 0.3 - 0.4, open models (models 4-9) are compatible with the data. This is 
our only robust conclusion about model viability. However, this is based on the assumption 
that there are no gross, unaccounted for, systematic uncertainties, and it only means that, 
in this case, the quoted error bars are not unreasonably small. (We emphasize that even if 

^ It might be signiflcant that for flat-A models the CMB anisotropy data seems to 

favour a larger h and a smaller O^, while some large-scale structure observations seem to 

favour a smaller h and a larger Qb (e-g-, Stompor et al. 1995; SSW). 

^ In our analysis here we have ignored upper limits. For the models we consider, the 

only seriously constraining upper limit is that of WDI (Tucker et al. 1993). This is mostly 

a serious constraint for the flat-A models (RS; RBGS), and is probably incompatible with 

these particular flat-A models (RS). 
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there are gross, unaccounted for, systematic uncertainties, the CMB anisotropy detections 
could stiU be compatible, but with different models.) 

It might be significant that a fairly large number of xfjj ^-re less than unity. While 
correlations between some data points certainly contribute to this, our understanding of 
the magnitude of the correlations leads us to suspect that mildly overconservative error 
bars on some of the data points might also be an issue. If this turns out to be more 
than just idle speculation, and if it can then be resolved, then, in combination with near- 
future improved small-scale CMB anisotropy data, and the tighter normalization error 
bars expected from the 4 year DMR data, our quantitative approach, based on using all 
available data, should allow for more robust conclusions about model viability. 

We are indebted to T. Banday, L. Page and J. Peebles, and also acknowledge useful 
discussions with E. Bertschinger, K. Gorski, G. Griffin, J. Gundersen, B. Netterfield, U. 
Seljak, and S. Tanaka. 
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TABLE 1 

Numerical Values for Model Parameters 



# 


Oo 


h 




(0)0 


0.1 


0.75 


0.0125 


(0)1 


0.2 


0.65 


0.0175 


(0)2 


0.2 


0.70 


0.0125 


(0)3 


0.2 


0.75 


0.0075 


(0)4 


0.3 


0.60 


0.0175 


(0)5 


0.3 


0.65 


0.0125 


(0)6 


0.3 


0.70 


0.0075 


(0)7 


0.4 


0.60 


0.0175 


(0)8 


0.4 


0.65 


0.0125 


(0)9 


0.4 


0.70 


0.0075 


(0)10 


0.5 


0.55 


0.0175 


(0)11 


0.5 


0.60 


0.0125 


(0)12 


0.5 


0.65 


0.0075 


(0)13 


1.0 


0.50 


0.0125 


(A)14 


0.1 


0.90 


0.0125 


(A)15 


0.2 


0.70 


0.0175 


(A)16 


0.2 


0.75 


0.0125 


(A)17 


0.2 


0.80 


0.0075 


(A)18 


0.3 


0.60 


0.0175 


(A)19 


0.3 


0.65 


0.0125 


(A)20 


0.3 


0.70 


0.0075 


(A)21 


0.4 


0.55 


0.0175 


(A)22 


0.4 


0.60 


0.0125 


(A)23 


0.4 


0.65 


0.0075 


(A)24 


0.5 


0.60 


0.0125 


(A)25 


1.0 


0.50 


0.0125 
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FIGURE CAPTIONS 

Fig. 1.- Five reduced (for 32 degrees of freedom) computed for the nominal value 

of the DMR normalization, and accounting for model "errors" , as a function of model 
number. The vertical short-dashed line divides the open models from the flat-A models. 
See Table 1 for model numbers and parameter values. The horizontal short-dashed lines 
are, in ascending order, at la (31.7% probability of xfjj being this large or larger), 2a 
(4.55% probability), 3a (2.70 x 10"^ probability), 4cr (6.33 x 10"^ probability), and 5a 
(5.73 X 10"'^ probability) for 32 de grees of freedom drawn from a gaussian distribution, 
and the horizontal long-dashed line is at the 95% probability level. 

Fig. 2.- As for Fig. 1, but computed for the nominal value of the DMR normalization, 

and now ignoring model "errors". Note that the vertical axis scale is different. 

Fig. 3.- As for Fig. 1, but now computed for the upper la value of the DMR normalization, 

and accounting for model "errors". Note that the vertical axis scale is different. 

Fig. 4.- As for Fig. 1, but now computed for the upper la value of the DMR normalization, 

and ignoring model "errors". Note that the vertical axis scale is different. 

Fig. 5.- As for Fig. 1, but now computed for the lower la value of the DMR normalization, 

and accounting for model "errors" . 

Fig. 6.- As for Fig. 1, but now computed for the lower la value of the DMR normalization, 
and ignoring model "errors" . 
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